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ABSTRACT 

One  of  the  most  striking  manifestations  of  instability  in  solid  mechanics  is  the  localiza¬ 
tion  of  shear  strain  into  narrow  bands  during  high  speed,  plastic  deformations  of  metals. 
According  to  one  theory,  the  formation  of  shear  bands  is  attributed  to  effective  strain¬ 
softening  response,  which  results  at  high  strain  rates  as  the  net  outcome  of  the  influence  of 
thermal  softening  on  the,  normally,  strain-hardening  response  of  metals.  Our  objective  is 
to  review  some  of  the  insight  obtained  by  applying  nonlinear  analysis  techniques  on  simple 
models  of  nonlinear  partial  differential  equations  simulating  this  scenario  for  instability. 
First,  we  take  up  a  simple  system,  intended  as  a  raradigm,  that  describes  isothermal  shear 
deformations  of  a  material  exhibiting  strain  softening  and  strain-rate  sensitivity.  As  it 
turns  out,  for  moderate  amounts  of  strain  softening  strain-rate  sensitivity  exerts  a  disspa- 
tive  effect  and  stabilizes  the  motion.  However,  once  a  threshold  is  exceeded,  the  response 
becomes  unstable  and  shear  strain  localization  occurs.  Next,  we  present  extensions  of  these 
results  to  situations  where  explicit  thermal  effects  are  taken  into  account. 
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One  of  the  most  striking  manifestations  of  instability  in  solid  mechanics  is  the  localiza¬ 
tion  of  shear  strain  into  narrow  bands  during  high  speed,  plastic  deformations  of  metals. 
According  to  one  theory,  the  formation  of  shear  bands  is  attributed  to  effective  strain¬ 
softening  response,  which  results  at  high  strain  rates  as  the  net  outcome  of  the  influence  of 
thermal  softening  on  the,  normally,  strain-hardening  response  of  metals.  Our  objective  is 
to  review  some  of  the  insight  obtained  by  applying  nonlinear  analysis  techniques  on  simple 
models  of  nonlinear  partial  differential  equations  simulating  this  scenario  for  instability. 
First,  we  take  up  a  simple  system,  intended  as  a  paradigm,  that  describes  isothermal  shear 
deformations  of  a  material  exhibiting  strain  softening  and  strain-rate  sensitivity.  As  it 
turns  out,  for  moderate  amounts  of  strain  softening  strain-rate  sensitivity  exerts  a  dissipa¬ 
tive  effect  and  stabilizes  the  motion.  However,  once  a  threshold  is  exceeded,  the  response 
becomes  unstable  and  shear  strain  localization  occurs.  Next,  we  present  e.xtensions  of  these 
results  to  situations  where  explicit  thermal  effects  are  taken  into  account. 


INTRODUCTION 

Shear  bands  are  regions  of  intensely  concentrated  shear  strain  that  are  observed  during 
the  plcistic  deformation  of  many  materials  and  often  precede  rupture.  Their  occurrence 
is  typically  associated  with  strain  softening  type  response,  past  a  critical  strain,  of  the 
measured  average  shear  stress  versus  the  measured  averaged  shear  strain.  The  diversity 
of  materials  and  deformation  conditions  associated  with  observations  of  shear  bands  has 
lead  to  various  theories  and  corresponding  models  for  explaining  their  formation.  Some 
type  of  strain-softening  mechanism  lies  in  the  core  of  most  theories,  however,  the  origin  of 
softening  is  often  context  dependent,  hinging  on  the  material  under  consideration  and  the 
geometric  and  loading  circumstances. 

In  this  article,  we  focus  on  a  theory  suggested  for  explaining  the  formation  of  shear 
bands  during  high  speed,  plastic  deformations  of  metals.  It  was  recognized  by  Zener 
and  Hollomon  [34]  that,  at  high  speed  processes,  the  effect  of  the  deformation  speed  is 
twofold:  First,  an  increase  in  the  deformation  speed  changes  the  deformation  conditions 
from  isothermal  to  nearly  adiabatic.  Second,  strain  rate  has  an  effect  per  se  and  needs  to 
be  included  in  the  constitutive  modeling. 

Under  isothermal  conditions,  metals  in  general  strain  harden  and  exhibit  a  stable 
response.  As  the  deformation  speed  increases,  the  heat  produced  by  the  plastic  work 
triggers  thermal  effects.  In  particular,  thermal-softening  properties  of  metals  may  outweigh 
the  tendency  of  the  material  to  harden,  so  that  the  combined  outcome  results  to  (effective) 
softening.  A  destabilizing  feedback  mechanism  is  then  induced,  operating  according  to  the 
following  scenario  (  Clifton,  Duffy,  Hartley  and  Shawki  [11]):  Nonuniformities  in  the  strain 
rate  result  in  nonuniform  heating.  Since  the  material  is  softer  at  the  hotter  spots  and 
harder  at  the  colder  spots,  if  heat  diffusion  is  too  weak  to  equalize  the  temperatures,  the  .  codes 
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initial  nonuniformities  in  the  strain  rate  are,  in  turn,  amplified.  This  mechanism  tends  to 
localize  the  total  deformation  into  narrow  regions.  On  the  other  hand,  there  is  opposition 
to  this  process  by  ’’viscous  effects”  induced  by  strain-rate  sensitivity.  The  two  effects  are 
competing  and  which  one  prevails  depends  on  the  relative  weights  of  thermal  softening, 
strain  hardening  and  strain-rate  sensitivity,  as  well  as  the  loading  circumstances. 

Experimental,  numerical  and  linearized  analysis  studies  indicate  that,  at  least  when 
the  degree  of  thermal  softening  is  large,  the  competition  results  to  instability  in  the  form  of 
strain  localization  and  formation  of  shear  bands  {e.g.  Clifton,  Duffy,  Hartley  and  Shawki 
[11],  Drew  and  Flaherty  [17],  Wright  and  Batra  [31],  Burns  [3],  Anand,  Kim  and  Shawki 
[1],  Wright  and  Walter  [32]).  There  is  an  extensive  literature  on  the  subject  and  the 
reader  is  referred  to  the  article  by  Shawki  and  Clifton  [23]  for  an  excellent  survey.  In  the 
present  article  we  review  recent  analyticcd  results  on  simple  constitutive  models,  where  the 
interplay  of  various  contributing  factors  on  shear  strain  localization  is  assessed. 

DESCRIPTION  OF  THE  MODEL 

Typically,  shear  bands  appear  and  propagate  as  one-dimensional  structures  (up  to  interac¬ 
tion  times).  Focusing  on  their  development  and  evolution,  we  consider  a  one-dimensional, 
simple  shearing  deformation.  In  a  Cartesian  coordinate  system,  consider  an  infinite  plate 
of  unit  thickness  located  between  the  planes  x  =  0  and  i  =  1.  For  a  simple  shearing 
motion,  the  only  nonvanishing  velocity  component  corresponds  to  the  shearing  direction 
and  all  the  field  variables  depend  only  on  the  x-coordinate.  Neglecting  the  normal  stresses, 
the  thermomechanical  process  is  described  by  the  list  of  variables: 

velocity:  v{x,t) 
shear  strain  :  u{x,t) 
temperature:  6{x,t) 
heat  flux  :  q{x,t) 
shear  stress  ;  <7(x,t). 

They  are  connected  through  the  balance  of  linear  momentum 

(1) 

kinematic  compatibility 

ut=v^,  (2) 

and  the  balance  of  energy  equations 

+  crut .  (3) 

Above,  the  reference  density  and  specific  heat  have  been  taken  constants  (set  equal  to  1). 
The  following  set  of  constitutive  assumptions  are  used  throughout: 

9  =  0,  (4) 
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a  =  f{d,u,ut),  (5) 

where  /  is  a  smooth  function  with  /(0,  u,0)  =  0  and  fp(B,u,p)  >  0,  for  p  ^  0.  A  discus¬ 
sion  of  the  constitutive  assumptions  follows,  but  note  that,  in  terms  of  classification,  the 
resulting  system  (1-5)  belongs  to  the  framework  of  one-dimensional  thermoviscoelasticity 
and  is  compatible  with  the  requirements  imposed  by  the  Clausius-Duhem  inequality. 

It  is  instructive  to  present  a  derivation  of  (4  —  5),  starting  from  a  constitutive  theory 
appropriate  for  thermal  elastic-viscoplastic  materials  (see  Shawki  and  Clifton  [23]  for  a 
detailed  derivation  starting  from  a  three-dimensional  model;  also  Needleman  [22]).  In  this 
context,  it  is  assumed  that  the  shear  strain  u  be  decomposed,  additively,  into  elastic  and 
pleistic  components:  u  =  +  u^.  The  elastic  component  u®  satisfies  linear  elasticity  with 

shecir  modulus  G,  that  is  u®  =  The  evolution  of  the  plastic  component  is  dictated  by 
a  plastic  fiow  rule: 

u^t  =  g{6,uP,(7)  or  a  =  f{9,u^ ,u^)  ,  (6) 

where  g  is  an  increasing  function  in  the  variable  <t,  and  f[9,  u,  •)  is  the  inverse  function  of 
g(9,  w,  •).  In  summary, 

-f  =  u 

1  ^  (7) 

—at  +  g(9,uP,a)  =  . 

For  the  heat  flux  either  a  Fourier  law,  q  =  k9x,  is  used,  or  it  is  asserted  that  the  process 
is  adiabatic,  i.e.  ^  =  0.  Imposing  adiabatic  conditions  projects  the  belief  that,  at  high 
strain  rates,  heat  diffusion  operates  at  a  slower  time  scale  than  the  one  required  for  the 
development  of  a  shear  band.  It  appears  a  plausible  assumption  for  the  shear  band  initi¬ 
ation  process,  but  not  necessarily  for  the  evolution  of  a  developed  band,  due  to  the  high 
temperature  differences  involved.  Equations  (1),  (3)  and  (7),  together  with  a  Fourier  law 
for  the  heat  flux,  are  regarded  as  modeling  the  essential  features  of  shear  strain  localization 
in  thermal  elastic- viscoplastic  materials.  They  lead  to  the  model  (1  —  5),  by  neglecting  the 
elastic  effects  {G  —  oo  and  u  =  u^)  and  working  under  adiabatic  conditions. 

Viewing  (5)  as  a  plastic  flow  rule  suggests  the  following  terminology:  The  ma'.erial 
exhibits  thermal  softening  at  the  state  variables  (0,  u,p)  where  fe{9,u,p)  <  0,  strain  hard¬ 
ening  where  fu{9,u,p)  >  0  ,  and  strain  softening  where  fu{9,u,p)  <  0.  Th*-  amounts 
of  the  slopes  of  /  in  the  directions  9,  u  and  p  measure  the  degree  of  therm  J  softening, 
strain  hardening  (or  softening)  and  strain-rate  sensitivity,  respectively.  The  difficulty  of 
performing  high  strain-rate  experiments  causes  uncertainty  as  to  the  sppcific  form  of  the 
constitutive  relations.  As  a  result  various  models  have  been  used  in  the  literature.  Ex¬ 
amples  of  specific,  experimentally  fitted,  constitutive  models  can  be  found  in  Shawki  and 
Clifton  [23].  The  following  are  used  for  the  purpose  of  illustration  in  this  article: 
Empirical  power  law  with  parameters  <  0,  and  n  >  0 


tf  XL  • 


Exponential  type  of  law 


a  =  e 
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The  pax2imeters  u,  k  and  n  in  the  power  law  (8)  serve  as  measures  of  the  degree  of  thermal 
softening,  strain  hardening  (or  softening)  and  strain-rate  sensitivity.  The  exponential  law 
(9)  does  not  exhibit  any  strain  hardening  and  the  parameters  0  and  n  measure  the  degree 
of  thermal  softening  and  strain-rate  sensitivity,  respectively. 

Isothermal  vs.  adiabatic  deformations 

To  illustrate  the  effect  of  thermal  softening  on  spatially  uniform  deformations,  the  isother¬ 
mal  and  adiabatic  cases  are  contrasted.  Consider  a  theoretical  experiment  where  the 
plate  is  subjected  to  steady  shearing.  Mathematically,  that  corresponds  to  prescribing  the 
boundary  velocities,  say  u  =  0  at  i  =  0  and  t;  =  1  at  x  =  1. 

(i)  In  an  isothermal  deformation  the  temperature  is  kept  constant,  say  6o,  by  appropriately 
removing  the  produced  heat  due  to  the  plastic  work.  The  “measured”  stress-strain  response 
in  this  ideadized  situation  coincides  with  the  T,—U  graph  of  S  =  /(^o ,  U,  1).  The  slope  of  the 
graph  is  measured  by  /u(^0i  U,  1);  hence,  for  a  stradn-hardening  material  E  is  monotonically 
increasing  with  U. 

(ii)  The  situation  in  an  adiabatic  deformation  is  understood  by  studying  a  special  class  of 
solutions  to  (1  —  5)  describing  uniform  shearing.  These  are 


where 


v(x,t)  =  X 

u(x,t)  =  U(t)  :=  t  +  uo 
6{x,t)  =  0(t), 


^  I) 


(10) 


(11) 

0(0)  = 

and  Uo,  6o  are  positive  constants,  standing  for  the  initial  strain  and  temperature.  Accord¬ 
ingly,  the  resulting  stress  is  given  by 


^(0  “  "b  “0i  1) 


(12) 


The  effective  stress-strain  curve  coincides  with  the  E  —  t  graph,  and  the  material  exhibits 
effective  hardening  in  the  increasing  parts  of  the  E  —  <  graph  and  effective  softening  in  the 
decreasing  parts.  Thus,  pending  on  the  sign  of  the  quantity  fe  f  +  fw,  the  combined  effect 
of  strain  hardening  and  thermal  softening  can  deliver  net  softening.  For  instance,  consider 
the  case  of  a  strain-hardening  {k  >0)  power  law  (8).  The  uniform  shearing  solutions  then 
read 

U {t)  =  t  Uq 


0(0  = 


0^0-''  + 


1  -  1/ 


k  + 

E(0  =  0‘'(O(<  +  uo)" 

and  a  simple  computation  yields 


(13) 


dt 


=  0(O""~'(^  +  uo)"" 


2  k 


u  +  k 


+ 


k 


LI--I-1  (t-buo)^+i 


0^-'' 

’FTT' 


(14) 
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For  parameter  values  reinging  in  the  region  v  +  k  <  0,  E(t)  may  initially  increase  but 
eventually  decreases  with  t. 

A  hierarchy  of  models 

It  is  generally  maintained  that  strain  softening  has  a  destabilizing  influence,  tending  to 
amplify  small  nonuniformities.  At  the  analytical  level,  this  belief  is  illustrated  by  the 
simple  model 

vt  =  r(u)i 

^  ’  (15) 

Ut-Vx, 

with  r'(u)  <  0,  describing  isothermal  motions  of  a  strain-softening,  inelastic  material. 
The  system  (15)  is  elliptic  in  the  t-direction,  and  the  initial  value  problem  is  ill-posed. 
Nevertheless,  it  admits  the  class  of  uniform  shearing  solutions  v{x,t)  =  x,  u(x,t)  =  t  -1-  uq. 

Because  of  the  inherent  instability  induced  by  strain  softening,  it  has  been  postulated 
that  higher  order  effects,  such  as  strain-rate  dependence,  are  triggered  in  the  course  of 
the  motion  and  need  to  be  accounted  for  in  the  constitutive  modeling  {e.g.  Needleman 
[22],  Wu  and  Freund  [33]).  Strain-rate  dependence  competes  with  strain  softening,  tending 
to  diffuse  nonuniformities  in  the  strain-rate  and/or  the  stress,  and  it  may  hinder  or  even 
altogether  suppress  instability  [26,  2S].  First,  we  present  a  quantitative  analysis  of  this 
competition  in  the  context  of  the  isothermal  model 

Ut  =  (r(u)Ox 


Ut  =  Vj 


with  the  function  t{u)  satisfying 


r(u)  >  0  and  r'(u)  <  0.  (Ho) 

Note  that  (16)  is  a  regularization  of  (15)  incorporating  rate-dependence  effects.  Different 
types  of  rate-dependent  constitutive  laws  have  also  been  employed  in  the  studj'  of  shear 
bands  (e.g.  Wu  and  Freund  [33]),  as  well  as  strain- gradient  dependent  constitutive  laws 
(e.g.  Coleman  and  Hodgdon  [12],  Zbib  and  Aifantis  [35]). 

Next,  we  proceed  in  the  framework  of  the  thermomechanical  model  (1  —5),  and  review 
investigations  [27,  30]  of  the  model 

Vt  =  {^{6yu)v^)x 

Ut  =  Vx  (17) 

Ot  =  fi{6,u)v"'^^ 

System  (17)  arises  by  substituting  in  (1  —  4)  the  simplified  version  of  the  constitutive 
relation  (5) 

a  =  fi{6,u)ut,  (18) 


where  the  smooth  function  fi(6,u)  sati.sfies 


fi(0,u)  >  0, 


(^i) 
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rendering  (17)i  parabolic,  and 


fie{e,u)<0,  (H2) 

determining  a  thermally  softening  material.  Both  strain  hardening  >  0),  strain 

softening  u)  <  0),  or  combinations  thereof  are  admitted  as  possibilities.  Instead,  it 

is  postulated  that 

6fj.eiB,u)  +  unu(6,u)  <0.  (H3) 

This  guarantees  that,  in  the  course  of  the  motion,  thermal  softening  prevails  over  strain 
hardening  so  that  the  material  eventually  exhibits  net  softening. 

When  is  independent  of  6,  the  first  and  second  equations  in  (17)  decouple 

from  the  third  and  lead  to  (16).  The  system  (16),  when  supplemented  with  Hypothesis 
(Ho),  may  be  viewed  as  an  isothermal  model  where  thermal  effects  are  implicitly  taken 
into  account.  That  is,  the  combined  outcome  of  thermal  softening  and  strain  hardening,  a 
coupling  effected  through  the  neglected  energy  equation,  results  to  net  softening,  which  is 
in  turn  reintroduced  as  an  assumption  by  imposing  (Hq). 

Another  interesting  simplification  occurs  when  fi{6,u)  is  independent  of  u,  corre¬ 
sponding  to  materials  that  do  not  exhibit  strain  hardening.  In  this  case  the  first  and  third 
equations  in  (17)  decouple  from  the  middle  one  and  lead  to  the  system 

V,  = 

e,  =  mo;*'  ■ 


Note  that  (19)  can  also  be  interpreted  as  describing  simple  shearing  of  a  non-Newtonian 
fluid  with  temperature  dependent  viscosity.  Strain  independent  constitutive  relations  were 
the  first  to  be  analytically  investigated  [15,  5-7,  25,  8,  2],  and  their  study  motivated  several 
of  the  ideas  presented  here. 

Initial  and  boundary  conditions 

The  governing  system  of  partial  differential  equations  is  taken  over  (x,<)  €  [0, 1]  x  {f  >  0). 
It  is  supplemented  with  initial  conditions 

v(x,0)  =  vo(x) ,  u(x,0)  =  uo(x),  B(x,0)  =  Bo(x) ,  (20) 

for  0  <  X  <  1,  and  as  a  consequence 

<7(x,0)  =cro(x)  := /z(0o(3:),uo(x))uo^(x).  (21) 

On  occasion,  the  initial  data  arc  used  to  provide  initial  nonuniformities  in  one  or  more  of 
the  field  variables  in  order  to  monitor  their  evolution. 

Two  distinct  sets  of  boundary  conditions  are  considered: 

a(0,t)  =  <r(l,0  =  l,  00,  (22)s 

corresponding  to  prescribed  tractions  at  the  boundaries,  and 

i;(0,f)  =  0,i;(l,f)  =  l,  f>0,  (22)v 
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corresponding  to  steady  shearing  of  the  plate  boundaries.  The  resulting  initial-boundary 
value  problems,  consisting  of  (17)  (or  (16),  or  (19)),  (21)  and  (22),  are  denoted  by  {V)s 
when  (22)5  are  used,  and  by  (V)v  when  (22)^  are  used. 

Preliminary  information 

For  the  validity  of  the  model,  it  is  necessary  that  solutions  comply  with  the  sign  restrictions 
(T>0,«t>0,«>0  and  0  >  0.  Indeed,  the  constitutive  law  (16)  arises  as  a  simplification 
of  (7)  by  neglecting  the  elastic  effects;  on  the  other  hand,  when  ut  becomes  zero  the  plastic 
flow  stops  and  elastic  effects  become  dominant.  The  last  two  restrictions  are  then  a  matter 
of  a  transformation  of  dependent  variables. 

An  existence  and  continuation  theory  that  covers  the  above  models  is  carried  out  in 
Ref  [28].  It  requires  smooth  initial  data,  together  with  hypothesis  {Hi ),  and  yields  a  unique 
classical  solution  (u(r,  t),  u(r,  t),  d(x,  <))  of  (P)s  (or  (P)v),  satisfying  (T>0,ut>0,u>0 
and  0  >  0.  The  solution  is  defined  on  a  maximal  time  interval  of  existence  [0,r*),  which 
in  general  could  be  infinite  or  finite.  An  existence  theorem  of  weak  solutions  for  (19)  (with 
n  =  1)  is  carried  out  by  Charalcimbakis  and  Murat  [8]. 

In  the  remainder  we  record  information  on  the  behavior  of  solutions,  with  the  following 
broad  objectives:  (a)  To  assess  in  a  quantitative  fashion  the  interplay  of  thermal  softening, 
strain  hardening  and  strain-rate  sensitivity  and  their  effect  on  the  response  of  shearing 
motions,  (b)  To  examine  the  circumstances  that  lead  to  the  developement  or  evolution  of 
nonuniformities  in  the  field  variables. 

ON  THE  COMPETITION  OF  STRAIN  SOFTENING  AND  STRAIN-RATE 
SENSITIVITY 

First,  we  take  up  the  isothermal  model  (16),  associated  with  the  constitutive  relation 

cr  =  r(u)  u”  (23) 

under  hypothesis  (i7o).The  system  (16)  may  be  viewed  for  small  n  as  a  regularization  of 
the  elliptic  system  (15).  On  the  other  hand,  (16)  belongs  to  what  is  formally  classified  as 
hyperbolic-parabolic  systems,  what  appears  in  discord  with  the  view  as  a  regularization 
of  an  elliptic  system.  To  reconcile  the  two  aspects,  recall  that  (16)  admits  the  uniform 
shearing  solutions 


v(x,t)  =  X  u(x,f)=t  +  uo  d-(x,t)  =  T(t  +  Uo)  , 

and  consider  a  perturbation  of  them 

v(x,t)  =  X  +  V(x,t) 
u(x,t)  =  (t  +  Uo)  +  U(x,t) . 

The  linearized  equations  for  (F,  U)  are  eeisily  computed: 

Vt  —  n  r(<  +  Uo)  Vri  -|-  T\t  uo)  U z 
Ut  =  Vz 


(24) 


(25) 


(26) 
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Under  hypothesis  (Hq),  when  n  =  0,  (26)  becomes  an  elliptic  system. 

Regarding  the  effect  of  the  boundary  conditions,  we  remark  that,  for  a  strain  softening 
material,  the  amount  of  external  work  required  to  prescribe  the  tractions  is  greater  than 
the  amount  required  to  prescribe  the  velocities.  It  is  thus  expected  that  the  shearing 
deformation  is  more  intense  in  the  former  case,  causing  favorable  conditions  for  unstable 
response.  In  the  course  of  shearing,  under  either  kind  of  loading,  the  material  is  strained 
and  the  diffusion  coefficient  in  the  parabolic  equation  (16)i  is  decreasing.  If  the  decrease 
is  too  rapid  and/or  spatially  nonuniform,  it  is  conceivable  that  the  diffusion  is  unable  to 
stabilize  the  motion,  possibly  resulting  to  spaticdly  nonimiform  structures. 

Both  of  the  initial-boundary  value  problems  {'P)s  and  {'P)v  for  the  quasilinear  system 
(16)  admit  classical  solutions  (i;(r,t),u(r,t)),  defined  on  a  maximal  time  interval  [0,T*) 
and  subject  to  the  sign  restrictions  (7>0,itt>0,u>0  [28].  As  an  outcome  of  a 
continuation  theorem,  if  T*  is  finite,  at  the  critical  time  the  strain  and  strain  rate  blow  up 

lim  sup  u(x,t)  =  oo 
lim  sup  sup  ut{x,t)  =  oo 

T*  0<r<l 

in  a  way  that  the  total  stress  remains  bounded 

0  <  <7  =  t(u)u?  <  max  ao(x).  (28) 

The  response  predicted  by  (27)  and  (28)  is  compatible  with  expectations  that  the 
strain  and  strain-rate  are  very  large  in  the  interior  of  a  shear  band  and  suggests  to  look  at 
shear  band  formation  as  a  blow-up  type  of  problem.  This  point  of  view  appears  instructive 
as  a  first  attempt  of  analytical  approach  to  the  problem,  but  is  probably  too  restrictive. 
It  leads  though  to  the  question  of  characterizing  whether  solutions  are  globally  defined  or 
break  down  in  finite  time. 

To  this  end,  it  is  useful  to  write  (23)  in  the  form 

$(u(x,<))  =  $(uo(a;))  +  /  crn(x,r)(ir,  (29) 

Jo 

where 

i{u):=  J'riO-di.  (30) 

Then,  blow  up  of  solutions  is  related  to  the  convergence  of  the  integral 

/oo 

(31) 

Indeed,  in  the  case  of  the  problem  ('P)s,  evaluating  (31)  at  x  =  0  or  x  =  1  leads  to 

$(u(f,f))  =  $(«o(0)+<  ,  2  =  0,1.  (32) 
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In  turn,  if  $((X))  <  oo,  (32)  implies 

u(i,t) oo  as  t-*Ti,  (33) 

with  Ti  =  $(oo)  — $(«o(2))  <  oo  for  i  =  0, 1,  and  the  solution  breaks  down  in  finite  time.  It 
turns  out  that  this  criterion  provides  a  complete  characterization  for  blow-up  of  solutions 
in  the  case  of  (P}s 

T*  —  oo  if  and  only  if  $(oo)  =  oo  .  (34) 

However,  in  the  case  of  the  problem  (P)v, 


if  $(oo)  =  oo  then  T*  =  oo ,  (35) 

but  the  converse  is  an  open  problem. 

The  role  of  the  integral  (31)  in  localization  problems  is  brought  forth  in  the  work  of 
Molinari  and  Clifton  [21].  They  introduce  the  concept  of  ’’Too'localization”,  and  use  an 
approach  developed  by  Hutchinson  and  Neale  [19]  for  the  uniaxial  tensii'n  of  a  bai,  to 
characterize  ice-localization  of  the  strain  with  the  convergence  of  the  integral  (31). 

In  what  follows,  we  record  some  information  pertaining  to  the  behavior  of  solutions. 
The  two  problems  are  treated  separately,  (P)s  in  this  section  and  (P)v  in  the  following. 
The  reader  is  referred  to  [28,  29]  for  details  of  the  derivations. 

Prescribed  tractions 

The  class  of  positive,  decreasing  constitutive  functions  can  be  decomposed  into  two 
categories,  depending  on  whether  <l&(oo)  is  finite  or  infinite.  Roughly  speaking,  the  dividing 
line  consists  of  functions  r(u)  that  decay  to  zero  like  the  power  n“”. 

In  case  r(n)  =  tq  a  constant,  (r(x,t)  is  a  positive  solution  of 

(36) 

subject  to  the  boundary  conditions  (22)5,  and  a{x,t)  1  uniformly  in  a:  ^  [0, 1]  as  t  — >  oo. 
We  ask  whether  this  behavior  persists  for  positive  and  decreasing  functions  t{u). 

Two  representative  classes  of  functions  t{u)  are  considered:  Class  (He)  consists  of 
functions  that  decay  to  a  positive  constant  t(oo)  at  a  rate  dominated  by  a  power,  i.e.,  for 
some  c  >  0  and  a  >  0 

r(u)  >  r(oo)  >  0  ,  0  <  — r'(u)  <  ^  ,  u>0.  (He) 

Class  (Hp)  consists  of  powers 

r(u)  =  —  .  (Hp 

According  to  (34),  T*  is  infinite  when  m  <  n,  but  finite  when  m  >  n. 

The  system  (16)  can  be  recast  into  an  equivalent  formulation  consisting  of  a  reaction- 
diffusion  equation  (with  variable  diffusion)  coupled  to  an  ordinary  differential  equation 


/  /si  1  t'(^)  I  i\: 
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ut  =  — TT  V'JOJ 

T{U)n 

We  ask  the  question:  Which  out  of  the  reaction  and  diffusion  terms  of  the  parabolic 
equation  (37)  dominates  in  the  course  of  the  motion?  The  technical  vehicle  for  answering 
this  question  is  by  classifying  the  steady  super  and  sub-solutions  of  (37).  The  coefficients 
axe  viewed  as  unknown  functions  with  the  only  available  information  that  Ut  >  0.  For 
functions  of  class  (He)  the  diffusion  term  always  dominates,  while  for  powers  (Hp)  the 
diffusion  term  dominates  for  parameter  values  0  <  m/n  <  1/2  and  the  reaction  term 
dominates  for  m/n  >  1/2. 

Dominant  diffusion 

Combining  this  analysis  with  parabolic  type  energy  estimates  it  turns  out  that,  for  func¬ 
tions  r{u)  of  class  {He)  and  for  powers  {Hp)  satisfying  0  <  m/n  <  1/2,  the  shear  stress 
a{x^t)  is  attracted  to  the  constant  state  cr  =  1,  while  u{xff)  behaves  asymptotically  as  a 
function  of  time  alone  [28]:  As  f  — >  oo 

a{x,t)  =  \  +  0{t-^),  (39) 


j‘u(x,t)  ft 

I  r(^) "  =  < -f  0(  /  s~^ds) 

Juo(z)  Ji 


uniformly  on  [0, 1],  with  0  <  0  <  1.  Precisely,  /d  =  o  >  0  for  {He)  and  0  =  for  {Hp). 

The  response  in  this  case  is  similar  to  the  response  of  (36). 

Weak  diffusion 

For  powers  with  exponent  m/n  >  1/2  (including  the  parameter  values  where  solutions 
blow  up)  the  following  holds  [29]:  If  S{x)  is  a  positive  function  satisfying 

uo{x)^^-^S{x)i>^S..{x) 

m  (41) 

S(i)  >  cro(x) , 

then  the  corresponding  solution  of  {V)s  satisfies 

cr{x,t)  <  S{x) .  (42) 

In  turn,  (29)  induces  corresponding  bounds  for  u{x,t).  The  implications  of  some  specific 
choices  for  S{x)  are  examined  below. 

Consider  first  the  one-parameter  family  of  functions  So(x)  =  1  —  Qrx(l  —  x),  with 
0  <  a  <  4.  Then  (41)  implies  the  restriction  on  the  initial  data 


arQ{x)  <  So{x)  ,  uo(x)  >  Ua  := 


n  2o 


,  0  <  X  <  1 . 
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Note  that  0  and  5a  — ♦  1  as  a  — >  0.  No  matter  how  close  to  the  state  cr  =  1  the 

initial  stress  is,  it  satisfies  thereafter 

cr(a:,<)  <  5a(ar)  =  1  -  a:i(l  -  x)  (44) 

and  can  no  longer  approach  the  state  cr  =  1.  We  conclude  that  in  the  region  m/n  <  1/2 
diffusion  is  too  weak  to  uniformize  the  stress,  and  initial  nommiformities,  however  small, 
persist  in  time. 

Next,  consider  a  positive,  convex,  and  symmetric  with  respect  to  x  =  |  function 
Sy(x),  that  drops  fast  from  5.^(0)  =  1  to  approximately  a  constant  level  7,  stays  there  over 
almost  the  entire  interval,  and  then  increases  fast  to  the  value  5-^(1)  =  1.  An  example  of 
such  a  function  is  provided  by 

5^(x)  =7  +  (l -7)[a;^  +  (l -x)^],  (45) 

for  N  a  large  positive  integer  and  0  <  7  <  1.  Let  the  initial  stress  uo(x)  be  symmetric 
with  respect  to  x  =  1,  with  Uo(x)  <  uo(0)  =  uo(l)-  Compatibility  with  (41)i  dictates 
that  uo(x)  be  sufficiently  large  in  a  neighborhood  of  x  =  0  and  x  =  1.  It  could  be  a 
sufficiently  large  constant.  Take  Sy(x)  as  the  initial  stress  and  uq(x)  as  the  initial  strain 
and  let  {(T^(x,t),u^(x,t))  be  the  corresponding  solution  of  {V)s- 

Then  satisfies  (42).  Regarding  the  behavior  of  u^,  we  look  at  two  separate 

regions: 

(i)  1/2  <  m/71  <  1:  Here  T*  =  00.  Also 


Uy(x,t)  < 


uj  "  (x) +(1  -  — )5-f  (x)< 
n 


n 

n  —  m 


(46) 


with  (46),  in  fact,  an  equality  at  x  =  0  and  x  =  1.  Since  uq  may  be  a  constant  function,  an 
initial  nonuniformity  in  the  stress,  no  matter  how  small,  may  induce  spatial  nonuniformities 
in  the  strain  u(x,t)  that  grow  in  time. 

(ii)  m/n  >  1:  Here  T*  <  00.  Also, 


u-f{x,t)  < 


n 


1  r  — /  — 1  ^  772  1 

5(n-y(x,t)  <  5-7  (x)  Ug  "  (x)  -  ( - l)5-7(x)t 

I  n 


(47) 


with  equalities  at  x  =  0  and  x  =  1.  By  the  choice  of  5-,(x)  and  ?<o(x),  the  right  hand  sides 
blow  up  for  the  first  time  at  the  boundary  points  x  =  0  and  x  =  1,  where  (47)  holds  as  an 
equality.  In  blowing  up,  they  appear  like  two  shear  bands  located  at  the  boundaries. 

The  above  considerations  suggest  for  the  case  of  prescribed  traction  loading:  First, 
that  it  is  possible  a  small  nonuniformity  in  the  stress  to  produce  large  strain  gradients 
in  either  long  time  (i),  or  even  shorter  times  (ii).  Second,  that  this  type  of  response  is 
concurrent  with  an  inability  of  the  material  to  diffuse  the  applied  stresses. 
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ON  THE  STABILITY  OF  THE  UNIFORM  SHEARING  SOLUTIONS 


In  this  Section,  we  take  up  the  initial-boundciry  value  problem  (P)k,  for  a  power  law  (ifp), 
and  investigate  the  stability  of  the  uniform  shearing  solutions  (24).  First,  we  present  a 
linearized  stability  analysis,  and  then  nonlinear  results  on  the  behavior  of  solutions. 

The  linearized  problem 

The  uniform  shearing  solutions  are  time-dependent,  thus  leading  to  linearized  equations 
with  variable  coefficients.  Investigations  that  account  for  the  time  dependence  of  the 
coefficients  in  the  linearized  equations  are  presented  by  Burns  [4],  Fressengeas  eind  Molinari 
[18]  and  Shawki  [24].  The  last  two  works  concern  models  that  include  explicit  thermal 
effects  through  the  energy  equation.  Fressengeas  and  Molinari  devise  a  relative  linear 
perturbation  analysis  that  examines  the  decay  (or  growth)  of  the  ratio  of  the  absolute 
perturbation  over  the  time-dependent  uniform  solution.  Shawki  proceeds  by  means  of  a 
quasistatic  approximation  for  the  linearized  equations  and  connects  the  issue  of  stability 
with  the  growth  in  the  kinetic  energy  of  the  perturbation  (see  article  in  this  volume). 

Here,  we  work  within  the  simple  context  of  (26)  for  a  power  law.  In  compensation, 
no  further  approximations  are  imposed  and  a  complete  picture  emerges.  It  is  convenient 
to  work  with  the  equations  for  the  linearized  displacement 

Ytt  =n(i  +  iio)“'"  Yizt  -  fn{t  +  uo)~'"“*  Yxx  ,  (48) 

with  Y(0,<)  =  Y(l,f)  =  0.  These  lead  to  (26)  by  setting 

V  =  Yt  17  =  Yx .  (49) 

The  general  solution  of  (48)  can  be  written  eis  a  linear  combination  of  eigensolutions 


Y{x,t)  =  Ck{t)  sin  knx  , 


fc=i 


where  Ck  are  Fourier  coefficients 


Ck{t)  =  2  I  Y{x,t)  sin  kTTx  dx 
Jo 


The  evolution  of  Ck  is  monitored  by  the  differential  equation 

(I’7r)2 


C\ck]  :=  Ck  + 


(nck - — — Ck)  =  0 . 

^  t  +  uo 


(t  -f  «o)' 


(50) 


(51) 


(32) 


For  simplicity,  suppose  that  y(a:,0)  =  ^>0(3^)  and  yi(x,0)  =  SVo{x),  where  6  >  0  is 
measuring  the  amplitude  of  the  initial  perturbation.  Then 


|cfc(0)l-h|cfc(0)l<6fc,  I- =  1,2,.... 


(53) 
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where  8k  =  0(8)  for  k  fixed,  and  its  dependence  with  respect  to  k  is  determined  by 
the  smoothness  of  Iq  and  Vq.  We  state  certain  facts,  that  are  proved  in  the  Appendix, 
concerning  the  response  of  solutions  of  (52)-(53)  for  0  <  m  <  1  and  n  >  0. 

(i)  Case  mjn  <  1. 

Given  any  e  >  0,  there  is  a  constant  Cc  independent  of  k  such  that 

\ck{t)\  <  C,  8k  (t  +  uq)^+‘  ,  \6k(t)\  <  C,  8k  (t  +  .  (54) 

The  above  estimate  is  sharp,  in  the  sense  that  there  are  initial  data  compatible  with  (53) 
such  that 

c;t(0  ^  (^  +  “o)”"  (55) 


(ii)  Case  m/n  >  1. 

For  any  initial  data  subject  to  (53),  there  is  a  positive  constant  C  independent  of  k  such 
that 

\ck{t)\  <  C  8k  (t  +  uo)  "  ,  |ct(<)|  <C8kit  +  uo)^-^  .  (56) 

Moreover,  there  are  initial  data,  compatible  with  (55),  such  that  for  any  £  >  0  there  is  a 
constant  Cg  with  the  property 


Ck(t)  >  Ct8k(t  +  uo)'^  '  .  (57) 

Since  the  constants  in  (54)  and  (56)  are  independent  of  k,  Y  and  V  =  Yt  inherit  these 
bounds,  with  the  same  time  rates  as  Ck  and  Ck  respectively.  Also  U  =  Yx  satisfies  a  bound 
with  the  same  time  rate  as  Y. 

To  summarize,  when  m/n  <  1  the  amplitude  of  perturbations  of  the  uniform  flow 
grows,  at  worse,  slower  than  the  uniform  flow.  By  contrast,  in  the  parameter  region 
m/n  >  1  perturbations  can  grow  faster  than  the  uniform  shear  flow,  and  the  uniform 
shearing  solutions  are  linearly  unstable. 

Nonlinear  response  for  prescribed  velocities 

The  linearized  stability  analysis  suggests  a  change  in  the  response  of  (V)v  across  the 
parameter  values  m/n  =  1.  This  is  confirmed  by  analysis  of  the  behavior  of  solutions  to 
the  nonlinear  problem  [28,  29]. 

Strong  strain-rate  dependence 

First,  we  present  a  nonlinear  stability  result  for  the  parameter  region  m/n  <  1.  Motivated 
by  the  form  of  the  uniform  shearing  solutions  and  the  scaling  invariance  of  the  underlying 
equations,  we  introduce  the  transformations 

v{x,t)  =  V{x,s{t)) 

u{x,t)  =  (t  +  l)U{x,s{t))  (58) 

a{x,t)  =  (t  +  1)“’”  E(x,s(<)) 
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with  s{t)  =  log(t  +  1).  A  straightforward  computation  shows  that  {T.{x,s),  U{x,  s))  satisfy 
the  system  of  reaction-diffusion  equations 

Consider  the  system  of  ordinary  differential  equations  obtained  by  neglecting  the 
diffusion  term  in  (59).  Their  integral  curves  are  the  lines  E  =  const.  U~’^  sketched  in 
Fig  1,  and  their  vector  field  vanishes  along  the  line  E  =  The  form  of  the  last  line 

changes  drastically  across  the  critical  parameter  value  mfn  =  1.  IfO<m/n<l,  the 
theory  of  Chueh,  Conley  and  Smoller  [10]  guarantees  that  (59)  admits  positively  invariant 
rectangles  of  arbitrary  size  in  the  first  quadrant  centered  around  the  line  E  =  (  see 

Fig  1).  (This  property  is  lost  in  the  complementary  region  m/n  >  1.) 


E 


—  m 


U 


(59) 


Fig.  1;  Invariant  regions  for  (59)  when  0  <  ^  <  1. 


Given  positive  initial  data,  let  U-,  [/+,  E_  and  E+  be  the  defining  coordinates  of  the 
smallest  rectangle  containing  the  curve  (<7o(x),  uo(i)),  0  <  x  <  1.  Then 

E_  <  E(x,s)  <  E+  ,  U-  <  U(x,s)  <  U+  .  (60) 

In  turn,  (58)  and  (60)  yield  the  relations; 

U-{t  -f  1)  <u(x,<)  <  U+{t  +  1) 

E-(f  -I- 1)"”  <<T{x,t)  <  E+(f  1)-'"  (61) 

Y.ZU2  <Ur(a',<)  <  . 
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Relations  (61)  provide  explicit  bounds  on  perturbations  of  the  uniform  shear  flow. 
Using  them  together  with  parabolic- type  energy  estimates,  it  turns  out  that,  for  m  < 
min{n,  1},  every  solution  is  asymptotically  attracted  to  a  uniform  shear  flow:  As  f  — >  oo. 


Vx{x,t)  =  1  + 

u(x,t)  =  t  +  0{t^)  (62) 

a(x,t)  =  f“"*(l  -f 

uniformly  on  [0,1],  with  P  =  max{^,m}  <  1.  The  role  of  the  constraint  m  <  1  is 
explained  in  [28].  We  remark  that  (62)  does  not  exclude  the  presence  of  nonuniformities 
in  the  strain  that  propagate  at  slower  time  rates,  and  is  thus  compatible  with  the  findings 
of  the  linearized  stability  analysis.  However,  as  indicated  by  the  bound  (61)3  on  the  strain 
rate,  no  catastrophic  growth  of  strain  is  predicted  in  this  parameter  range. 

Asymptotic  behavior  results  such  as  (62)  have  also  been  established  for  the  case  of  the 
constitutive  relation  cr  =  r(u)u(,  with  r(u)  increasing  up  to  a  critical  strain  and  decreaising 
(but  not  very  fast)  thereafter  [26].  However,  the  explicit  bounds  (61)  are  available  only  for 
a  power  law.  Such  results  are  interpreted  as  indicating  that  eventually  rate  dependence 
exerts  a  stabilizing  influence  that  washes  out  nonuniformities.  However,  at  intermediate 
stages,  there  may  be  local  amplification  of  nonuniformities  in  strain  connected,  for  instance, 
with  undergoing  through  a  sudden  drop  of  the  curve  t(u)  —  u. 

Weak  strain-rate  dependence 

In  the  region  m/n  >  1,  the  behavior  of  solutions  to  {V)v  is  quite  subtle  and  is  not  fully 
understood  at  the  present  time.  The  situation  in  the  case  of  prescribed  tractions  suggests 
that  unstable  response  be  expected,  at  least  for  large  values  of  the  ratio  m/n.  It  follows 
by  (2)  and  (22)  v  that  the  average  strain  grows  like  t 


(63) 


On  the  other  hand,  the  linearized  analysis  indicates  that  perturbations  of  the  strain  tend 
to  grow  faster,  like  t^.  It  is  conceivable  that  soon  large  nonuniformities  in  the  strain 
develop.  At  that  point  the  linearized  analysis  ceases  to  be  valid. 

We  discuss  now  the  effect  of  a  large  nonuniformity  in  the  strain  (see  [29],  also  Bertsch, 
Peletier  and  Verduyn  Lunel  [2]  for  a  similar  result  for  (19)).  The  analysis  applies  to 
constitutive  functions  that  include  powers  with  m/n  >  1  and  n  <  1.  Consider  the  initial 
data 

Vo(x)  =  X  0  <  X  <  1 

(  \  —  f  ^  0<a:<j/  —  5  ,  y  +  ^<3:<l  (64) 

~  \u(x)  y-8  <x  <y-\-6 

corresponding  to  a  concentration  of  strain  in  an  interval  of  length  28  around  a  fixed  point 
y,  and  a  linear  initial  velocity  profile.  We  are  thinking  of  tj(x)  as  a  bell  shaped  function 
that  attains  its  maximum  at  y,  say  U(y)  =:  Um,  and  satisfies  U(y  —  8)  =  U{y  +  8)  =  u 
and  V(x)  >  u  for  x  €  (y  —  ^ ,  y  +  ^)  (Fig  2). 
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It  turns  out  that  for  U m  sufficiently  large; 

(i)  either  T*  is  finite,  in  which  case  (27)  holds, 

(ii)  or,  as  f  — >  oo. 


X  6  [0,y  -  ^] 
X  €  [y  +  6, 1] 


(65) 


and  u{x,t)  converges  monotonically,  for  x  outside  the  band  (y-6,y  +  S),  to  some  bounded 
limiting  function. 


Fig.  2;  The  initial  data  uo{x),  vo(x)  in  (64). 


The  large  time  response  predicted  in  (65)  looks  like  a  fully  developped  shear  band,  where 
the  parts  of  the  material  to  the  left  and  right  of  the  band  have  unloaded  and  move  inde¬ 
pendently.  The  result  indicates  a  collapse  of  momentum  transfer  across  the  band. 

EXPLICIT  THERMAL  EFFECTS 

In  this  Section  we  present  information  regarding  the  behavior  of  solutions  of  (P)s  and  (V)v 
for  the  model  (17).  The  model  (17)  is  based  on  the  constitutive  relation  (18),  appropriate 
for  a  material  exhibiting  thermal  softening,  strain  hardening  and  strain-rate  sensitivity, 
and  is  studied  in  this  generality  in  Refs  [27]  (for  n  =  1)  and  [30]. 

The  analysis  of  the  uniform  shearing  solutions  (10)  indicates  that,  when  the  degree  of 
thermal  softening  is  large,  the  combined  effect  of  thermal  softening  and  strain  hardening 
results  to  net  softening.  Capturing  this  behavior  for  arbitrary  solutions  is  the  main  source 
of  analytical  difficulty  in  dealing  with  (17),  and  establishing  a  precise  analogy  among  the 
models  (17)  and  (16).  One  approach  to  accomplish  that  is  outlined  below.  Differentiating 
<rn  with  respect  to  t  one  obtains 

(<r")t  =  fi^(S,u)utt  -{-  ^fin-^(e,u)(dtfi(9,u))ut .  (66) 
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The  terrr  dtyi{d,u)  is  positive  or  negative  depending  on  whether  the  material  exhibits 
effective  hardening  or  softening  along  the  deformation.  Use  of  (66)  and  (17),  leads  to  an 
alternative  formulation  of  the  problem,  consisting  of 


(cr")t  =  /in(0,u)(Tij.  + 


n 


/i  +  i- 


{d,u) 


{fiei6,u)(T  +  fiui6,  u))  (cr  i  f 


(67) 


dt  =  (TUt  , 


(68) 


and 


<rn 


Ut  = 


/in(0,u) 


>  0, 


(69) 


As  in  the  case  of  (37-38),  we  ask  which  out  of  the  reaction  or  diffusion  terms  in  (67) 
dominates.  The  analysis  proceeds  by  classifying  steady  super  and  sub-solutions  of  (67). 
The  coefficients  are  viewed  as  unknown  functions  connected  through  (68),  which  controls 
the  relative  weight  of  temperature  and  strain  (see  [30]  for  details). 

This  approach  yields  a  characterization  of  the  behavior  of  the  field  variables  at  the 
critical  time,  when  T*  is  finite.  Hypotheses  (Hi  —  H3)  for  //(^, u)  are  used;  at  places  the 
additional  hypothesis 


/i”  =  00  forany/9>0,  {H4) 

is  employed.  If  T*  is  finite,  then  as  t  —*  T*  the  stress  remains  bounded  by  the  quantity 
S  =  max{maxo<r<i  o‘o{x) ,  maxo<r<i(^o((r)/wo((r))},  the  strain  and  strain  rate  behave  as 
in  (27),  while,  in  case  {H4)  also  holds,  the  temperature  satisfies 

lim  sup  ^(x,t)  =  oo.  (70) 

0<r<l 


Global  existence  or  blow-up  of  solutions  to  {V)s  or  {V)v,  is  related  to  the  convergence 
or  divergence  of  the  integreil 


(71) 


which  measures  how  fast  the  constitutive  function  pi  decays  (c/  {H3))  when  the  temper¬ 
ature  is  proportional  to  the  strain.  As  it  turns  out: 

(i)  For  either  of  the  problems  ('P)s  or  (V)v,  if  M{L)  =  00  then  T*  =  00. 

(ii)  For  the  problem  (7^)5,  if  M{k)  <  00  then  T*  <  00. 

Here  k  :=  min{l ,  ^o(0)/uo(0) ,  ^o(l)/t^o(l)}-  This  result  provides  a  characterization  of 
blow  up  for  stress  boundary  conditions,  but  only  sufficient  conditions  for  global  existence 
in  the  cjise  of  velocity  boundary  conditions. 

In  general,  the  above  criteria  depend  on  the  initial  data.  However,  for  certain  popidar 
models  used  in  the  study  of  shear  bands,  they  are  actually  independent  of  the  data. 
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(a)  For  a  power  law  (8), 


~n 

u+k+n 


<  OO 


ioT  u  +  k  >  —n 
for  u  +  k  <  —n 


(72) 


Hypothesis  {H^)  dictates  the  restriction  u  +  k  <  0,  while  hypothesis  (Hi)  is  fulfilled 
whenever  k  >  —n.  We  conclude  that  for  a  power  law  solutions  of  {V)s  and  {'P)v  exist  for 
all  times  for  parameter  values  —n  <  u  +  k  <  Q.  Solutions  of  ('P)s  blow  up  in  finite  time 
for  parameter  values  i/  +  fc  <  — n. 

(b)  For  an  exponential  law  (9)  with  ^  positive, 


M.{q)  =  — -  <  OO 

ap 


and  solutions  of  ('P)s  blow  up  in  finite  time  for  any  values  of  the  parameters. 

Power  law  and  prescribed  tractions 

In  the  case  of  a  power  law  and  for  shearing  deformations  caused  by  prescribed  tractions  it  is 
possible  to  obtain  a  better  understanding  of  the  behavior  of  solutions  [30].  The  parameter 
region  i/  -f  /:  <  0  can  be  decomposed  into  three  subregions  across  which  the  material 
response  changes  drastically.  The  situation  generalizes  the  results  for  strain  softening 
materials  (23)  of  class  (Hp). 

(i)  In  the  region  0  <  (i^  +  k)/n  <1/2  the  shear  stress  a{x,t)  is  attracted  to  the  constant 
state  cr  =  1,  as  t  — >  oo,  while  u(i,t)  and  9{x^t)  behave  asymptotically  as  functions  of  time. 

(ii)  In  the  region  1/2  <  (i/  +  A,’)/n  <  1  the  constant  state  <7  =  1  loses  its  stability  and 
spatial  nonuniformities  in  the  strain  can  develop  and  persist  in  time. 

(iii)  Finally,  in  the  region  {v  +  k)/n  >  1,  u{x,t)  becomes  infinite  in  finite  time. 

Most  interesting,  when  {u  +  k)fn  >  1/2,  there  are  initial  data  for  which  the  strain  u{x,t) 
develops  nonuniformities  around  x  =  0  and  x  =  1  and  looks  like  two  shear  bands  located 
at  the  boundaries.  This  response  is  concurrent  with  an  inability  for  diffusion  of  the  applied 
stresses. 

Strain  independent  materials 

At  the  present  time,  the  analysis  of  nonlinear  stability  of  the  uniform  shearing  solutions 
for  (17)  is  incomplete.  The  available  results  concern  the  special  cases  that  either  thermal 
softening  or  strain  dependence  are  absent.  For  a  strain  independent  power  law  (17)  takes 
the  form 

Vt  =  {e‘'v:h 

L,  (<4) 

01  =  ^  ^ 

with  u  <  0.  From  a  viewpoint  of  analysis  there  are  certain  similarities  between  (74)  and 
(16)  —  (Hp).  As  a  consequence,  the  response  of  solutions  to  (74)  parallels  the  results 
outlined  in  the  previous  section  for  (16). 

First,  if  1/  >  — n,  every  solution  is  asymptotically  attracted  as  t  >  oo  to  a  uni¬ 
form  shearing  solution  (Dafermos  and  Hsiao  [15],  Tzavaras  [25]).  Similarly  to  the  case  of 
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a  strain-softening,  rate-dependent  power  law,  a  derived  system  admits  invariant  regions 
(Bertsch,  Peletier  and  Verduyn-Limel  [2]  for  n  =  1,  Tzavaras  [30]  for  n  ^  1).  As  a  result, 
perturbations  in  the  temperature  and  strain  rate  of  the  uniform  flow  are  controlled  by 
bounds  (of  the  form  (61))  depending  on  the  data  and  evolving  according  to  the  time  rates 
of  the  uniform  flow,  given  in  (13)  for  =  0. 

Complementing  the  above,  Bertsh,  Peletier  and  Verduyn-Lunel  [2]  consider  the  case 
of  a  Newtonian  fluid  (n  =  1)  and  show  that,  for  parameter  values  v  <  —\  and  for  an  initial 
temperature  concentrated  in  a  narrow  region,  either  the  solution  of  (74)  blows  up  or  the 
velocity  has  an  asymptotic  profile  of  the  type  described  in  (65).  They  also  establish  the 
eisymptotic  behavior  of  solutions  in  the  borderline  case  u  — 

The  effect  of  other  types  of  loading  in  the  stability  of  shear  flows  for  temperatures 
dependent  fluids  is  studied  by  Charalambakis  [5-7|.  He  considers  loading  effected  by  mixed 
velocity-stress  boundary  shearing,  inertial  forces  or  periodic  boundary  shearing,  and  es¬ 
tablishes  sufficient  criteria  for  stability  of  such  flows. 

The  above  works  proceed  under  the  constitutive  assumption  (4)  of  a  non-conducting 
material.  Experience  accumulated  through  studies  of  initial  value  problems  in  thermome¬ 
chanics  (Dafermos  [13,  14],  Dafermos  and  Hsiao  [16])  suggests  that  heat  conduction  exerts 
a  dissipative  effect  hindering  unstable  response.  Due  to  the  high  temperature  differences 
across  a  shear  band,  heat  conduction  is  expected  to  play  a  significant  role  in  problems 
concerning  propagation  of  shear  bands,  and  may  conceivably  provide  shear  bands  with 
internal  structure. 

Insight  in  that  direction  is  offered  by  recent  works  (Chen,  Douglas  and  Malek-Madani 
[9],  Maddocks  and  Malek-Madani  [20]).  Their  analysis  applies  to  the  system  of  equations 
consisting  of  (1),  (3)  and  the  constitutive  relations 

q  =  kOj:  ,  a  —  e~^^Vi  (75) 

under  isothermal  boundary  conditions  for  the  temperature,  and  various  possibilities  of 
loading,  including  prescribed  tractions  and  prescribed  velocities.  The  existence  of  steady 
state  solutions  is  studied  and  the  presence  of  bifurcating  branches  of  solutions  is  established. 
In  addition,  the  stability  of  these  steady  states  is  studied,  using  a  direct  linearized  analysis 
and  energy  estimates  [9],  or  a  variational  stability  approach  [20]. 


APPENDIX 

We  present  here  the  technical  steps  that  provide  the  stability  analysis  for  the  linearized 
problem  (48).  It  is  based  on  properties  of  the  differential  operator 


C[g]  ■■=g  + 


(<  +  «o)’" 


m 

t  +  Uo 


9) 


with  fc  =  1,  2,  ... .  Note  that  if  m/n  =  1  then  C[{t  +  uq)]  =  0. 

Lemma.  Suppose  that  0  <  m  <  1  and  n  >  0. 

For  parameter  values  m/n  >  1: 

(ii)  If  /+(t)  =  {t  +  uo)^  then  C[f+]  >  0. 
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(i2)  Given  any  £  >  0,  there  is  a  positive  constant  A'j  depending  on  £  but  not  on  k,  such 
that  the  function 

/-(t)  =  (<  +  uo)^-'  +  A. 

satisfies  £[/-]  <  0,  for  <  >  0. 

For  parameter  values  mjn  <  1: 

(111)  There  is  0  <  q  <  1  such  that  the  function  H+(t)  =  (t  +  iio)“  satisfies  C[H+]  >  0,  for 
t  >  0. 

(112)  Given  any  £  >  0,  there  is  a  positive  such  that  the  function 

MO  =  (f  +  uo) " 


satisfies  >C[/i-(-]  >  0,  for  t  >Tc. 

(iia)  If  h-{t)  =  {t  +  uq)^  then  <  0. 

Proof .  First,  observe  that 


—  ^  ^  t  f  if  _  -v  ^  V  O 

P[{t  +  Uq)  n  ]  =  - ( - l)(t  +  Uo)  " 

n  n 

is  positive  for  m/n  >  1  and  negative  for  m/n  <  1.  This  shows  (ii)  and  (iis). 

Next  we  turn  to  (i2).  A  straightforward  but  lengthy  computation  yields 

£[/_  I  =  (f  +  ^[A  -  en(A-:r)^((  +  uo)‘-”>]  -  h\  ■ 

where  A  =  ( ^  —  £)(  ^  —  £  —  1).  If  A  is  negative  the  result  follows.  So,  we  consider  the  case 
0  <  £  <  ^  —  1,  for  which  A  is  strictly  positive.  Let  Tj  be  such  that  {Tg  +  uo)^“'"  >  Xfen-'^. 
Then  £[/-]  <  0  for  t  >  Tg.  For  the  remaining  interval  [0,T’j]  we  attain  C[f-]  <  0,  by 
choosing  the  positive  constant  Kg  sufficiently  large. 

Finally  we  turn  to  (iii)  and  (ii2)-  A  computation  yields 

C[  (t  +  uo)“  ]  =  (^  +  [(no  -  m){kw)'^  {t  +  -  o(l  -  «)]  • 

Let  first  o  be  sufficiently  close  to  1,  so  that  that  q(1  —  q)  <  7r^uo*~'”(na  —  m)  and 
m/n  <  o  <  1.  For  this  choice  of  o,  the  corresponding  function  H+{t)  satisfies  >  0 

for  t  >  0.  This  shows  (iii  )•  Set  now  o  =  m/n  +  e  and  let  Tg  be  such  that  (Tg  +  > 

(^  +  £)(1  —  ^  —  £)/£n7r^.  Then  for  t  >  Tg  the  function  h+{t)  satisfies  £[/r+]  >  0.  This 
shows  (ii2)  and  completes  the  proof  of  the  lemma. 

The  second  ingredient  of  the  linearized  analysis  is  a  comparison  estimate  for  the 
operator  C[g\.  Namely: 

(iii)  If  g{T)  >  0,  g(r)  >  0  and  C[g\  >  0  for  t  >  r,  then  g{t)  >  0  for  t  >  r. 

Indeed,  if  the  conclusion  is  violated  then  g  admits  a  positive,  local  maximum  at  some  point 
s  >  T.  At  s  it  is  g{s)  >  0,  g{s)  =  0,  ^(5)  <  0  and  thus  £[g](s)  <  0,  which  contradicts  the 
assumption. 
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Behavior  of  solutions  to  (52) 

Fix  k  a  positive  integer  eind  let  Ck(t)  be  a  solution  of  (52)  with  initial  data  subject  to 
(53).  Our  goal  is  to  monitor  the  time  evolution  of  Ck(t)  in  terms  of  the  size  Sk  of  the 
data.  In  the  sequel  the  ^-dependence  of  Sk  is  suppressed,  and  C  will  stand  for  a  generic 
positive  constant  that  is  independent  oi  t,  6  =  Sk  and  the  parameter  k.  Within  the  range 
0  <  m  <  1  and  n  >  0  two  separate  regions  are  considered. 

First  case:  mfn  <1 

Our  first  task  is  to  show; 

(iv)  If  Ck  satisfies  Cfc(O)  <  S,  Ck(0)  <  S  and  £[cjk]  =  0,  then  for  any  £  >  0  there  is  a  constant 
Cff  such  that 

Ckit)  <  C,S{i  +  ,  Cfc(t)  <Ce^(^  +  uo)^+-'■^  fc  =  l,  2,  ....  (a.l) 

To  this  end,  let 

g{t)  =  CiSH^(t)  -  Ckit) , 

where  H+(t)  is  as  in  (iii),  and  the  constant  Ci  is  chosen  sufficiently  large  so  that  <7(0)  >  0 
and  g(0)  >  J.  Since  C[g\  >  0  for  t  >  0,  (iii)  yields 

cjt(0  <  C\SHjf.{t)  =  Ci6(f  -f  uo)“  f  >  0 . 

To  improve  this  bound,  observe  that  (52)  implies 

{krtYn  .  {kTr)^m 

(t-l-uo)'”^*'  “  (t  -f 

SjkTTym 

~  ^\i  +  uo)^ 

Whenever  cjt  is  larger  than  CiSm/n,  Ck  is  decreasing.  Hence, 

Cifc(t)  <  max{cfc(0),  CiSm/n]  <  C2S  f  >  0. 

Consider  now  the  comparison  function 

g{t)  =  C3Sh^{t)  -  Ck{t) . 

In  view  of  (ii2),  there  is  Tg  so  that  C[g]  >  0  for  f  >  Tg.  Choose  Cj^g  =  Czis)  such  that 
g(Tg)  >  0  and  g(Tg)  >  0.  This  is  possible  because  of  the  preliminary  estimates  that  are  at 
our  disposal.  Then  (iii)  implies 

Ck{t)<C3,gS{t  +  Uo)^^‘  t>Tg.  (a.4) 

Combining  (a.2)  and  (a.4)  yields  the  first  part  of  (a.l). 


(a.2) 


(a.3) 
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To  show  the  second  part,  note  that  (52)  implies 


(k-rrYn  .  ^  „  c  {ki^Ym  ,  ,, 

Integrating  the  differential  inequality  and  performing  an  integration  by  parts  gives 

Ck{t)  <  Cfc(0)exp{  -Cik'^[{t  +  -  wo^“'”]  ) 

+  Cs.e  8  l\s  +  uo)  exp{  C^k^iis  +  uo)'"'"  -  (/  +  uo)'""*]  }  ds 

Jo 

<(5exp{-C4[(<  +  uo)'-"*  -uo'-"*]} 

C5,..5 


+ 


^  +  £-1 


n  iv“  ■  “u,  -u  exp{-C4^-^[(t +  uo)^  -  uq^  '"] } 

(1  -  m)C4  L 

+  (!---£)  f\s  +  uo)^+"-"exp{  C4t2[(5  +  uoY-^  -  (t  +  uoY-^]  }  ds 

Jo  1 

The  last  integral  is  estimated,  with  the  help  of  L’Hopital’s  rule,  to  arrive  at 

a(t)<C6,,6(<  +  uo) 


with  Cs.e  independent  of  k. 

Now  (54)  follows  by  applying  (iv)  consecutively  to  Ck  and  —Ck-  To  show  (55),  fix  the 
initial  data  to  be  c<,(0)  =  Ck(0)  =  S/2  and  test  the  comparison  function 

g{t)  =  Cfc(t)  -  C7Sh-{t) 

Then  C[g\  >  0  and,  upon  choosing  C7  appropriately  small,  we  will  have  $((0)  >  0  and 
g(0)  >  0.  It  follows  that  </(t)  >  0,  which  yields  (55). 

Second  case:  m/n  >  1 

A  similar  in  spirit  argument  zis  in  the  previous  case,  using  the  functions  f+{t)  and  /-(t), 
provides  (56)  and  (57). 
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